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Abstract 



We study the rare radiative dileptonic decays B°(B S ) — ► 7^ + .£ = e,/i) in the 
standard model. By using the B meson wave function constrained by non-leptonic 



decays, the branching ratios turn out to be of the order of 10 for B s — * 7/v/i~, 
7e + e~ and 10 -10 for B° —* 7 / u + /i~~, 7e + e~. Based on the study, these decays are 
accessible at the near future LHC-b experiment, which are useful to determine the 
B(B S ) wave function. 



1 Introduction 

The radiative and leptonic B decays have been the subject of many theoretical studies in the 
framework of the Standard Model (SM) and the search of new physics. These processes play 
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an important role in determining the parameters of the SM and some hadronic parameters in 
QCD, such as the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements, the meson decay 
constant fs, and fs s - Especially, the radiative leptonic decays can provide information on 
heavy meson wave functions [1] . Being rare decays in SM, they are also very sensitive to any 
new physics contributions. 

Due to the GIM mechanism [2 J , there is absence of flavor changing neutral current tran- 
sition at the tree level, thus, the pure leptonic processes B°(B S ) — > £ + £~ (£ = e, /i) can only 
occur via penguin and box diagrams. In addition, rare decays of heavy pseudoscalar meson 
into light lepton pairs are helicity suppressed. Some of the branching ratios are very small, 
for instance p] B r (B s — > = 2 x 1CT 9 and B r (B s — > e + e~) = 4 x 10~ 14 , so it is difficult 

for us to extract the B s meson decay constant information from these processes. For B° 
meson, the situation is even worse due to the smaller CKM matrix elements. Although the 
process B s — ► t + t~ is free from this helicity suppression mechanism, the branching ratio 
turns out to be around 8 x 10~ 7 in the SM [4J, which is much larger than the branching ratios 
of B s — > yU + /i~, B s — > e + e~, it is still hard to be detected at B-factory where the efficiency 
is not better than 10~ 2 . 

Because of massless neutrino, the decay B°(B S ) — > vv is helicity forbidden. While, if 
an additional real photon is emitted, the forbidden situation will change [5j [6]. Similarly 
the helicity suppression in the pure leptonic B decays B°(B S ) — > £ + £~ will be cured in 
the radiative decay B°(B S ) — > / y£ + £~~ . This is already shown in the simple constituent quark 
model calculation [7j, light cone QCD [8]. In this paper, we employ the B meson distribution 
amplitude derived from non-leptonic B decays to analyze these processes again. 

We organize our paper as follows: In Sec. 2, the relevant effective Hamiltonian will be given 
in the SM. In Sec. 3, the wave function of B meson will be used to calculate these processes, 
and later some comparison will be given. Finally, Sec. 4 includes a brief conclusion. 

2 Effective Hamiltonian 

Let us start with the quark level processes b q£ + £~, with q = s or d quark, £ = e or \i. 
The leading order Feynman diagrams are shown in FigHJ It is easy to see that the magnetic 



2 



penguin, Z penguin and box diagrams contribute to these processes. They are subject to 
QCD corrections which can be obtained by connecting the quark lines by gluon lines. The 
effective Hamiltonian in SM for them is [9]: 



n 



eff 



aGp 
V2n 



v tb v; q 



qP R (ff ~ l* 1 fib + C e 9 ff qYP L b 



where Pl = (1 — 75)/2, Pr = (1 +7 5 )/2, q 2 = (P + + P_) 2 , with P + and P_ are the momenta 
of lepton pair. C 7 , Cg , and C w are the QCD corrected Wilson coefficients, whose specific 
forms are given in Ref . [TU] . 
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Figure 1: Leading-order diagrams for b 



transition in SM. 



If an additional photon line is attached to any of the charged lines in Fig{H we will 
have the radiative leptonic decays b — > q^i^i' . In fact, there are two kinds of diagrams: 
photon connecting to the internal line of Fig.l, and photon connecting to the external line 
of Fig.l. For the first kind of diagrams, because of the effective operators are dimension-8 
instead of dimension-6, there will be a suppression factor of mj/m^, in the Wilson coefficients 
compared with the ones for b — > q£ + £~ . Therefore we will only consider the second category 
of diagrams. In this case, we only have dimension 6 operators with an additional photon 
from any of the fermion lines, which are shown in FigfSJ 

First, let us study the FigEte) and (f) with photon attached to external lepton lines. 
Being a pseudoscalar meson, B meson can only decay through axial current, so 

< 0\qa^P R b\B >= 0. (2) 
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Figure 2: Feynman diagrams for B q — > ^i^i in the SM (black dots are dimension 6 opera- 
tors). 



That is to say, in diagrams (e) and (f), the magnetic penguin operator O7 's contribution 
vanishes. Our numerical calculation shows that: the contribution from other operators for 
these two diagrams are also very small, we can neglect their contribution safely. Similarly for 
Figj2] (c) and (d), We find numerically they are also negligibly small comparing with Figj2] 

(a) and (b). Thus the main contribution to B°(B S ) — > r y^i~ comes from the Figj2] (a) and 

(b) , which agrees with the constituent quark model calculation [7]. 

Let us analyze the diagrams (a) and (b) in FigfSJ with photon emitted from the external 
quark lines h or q. After calculation, we get the amplitude: 



A 



2" 



e C 7 m b _ 
2 q l 



p 2 .k 
Pr 



*-fP L + >fP L . 



Pi-k 
}6 + m b 



b. 



C e 9 ff F ltl l + + C 10 F W 1 

f 2 + m 



(3) 



- Y 



Pi-k p 2 .k 
where pi, p 2 and k are the momenta of b, q quark and photon, respectively 
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3 Model calculations 



To simplify the decay amplitude in eq.([3]), we have to utilize the B meson wave function, 
which is not known from the first principal. Fortunately, many studies on non-leptonic B 
[TTj [T2"] and B s decays [T3] have constrained the wave function strictly: 



®B = ~Ei{i>B + m B )j 5 (j) B (x) 
V6 



where the distribution amplitude (fi B can be expressed as 

4>b{x) = N B x 2 {\ — x) 2 exp 



M 2 x 2 



2u 2 



It satisfies the normalization relation: 

• i 



B 



d> B (x)dx = , 



(4) 



(5) 



(6) 



with f B is the B meson decay constant. 

From definition of the wave function, we have pb = (1 — x)P B , p q — xP B , and the decay 
amplitude (J3J) is then: 



A = C 



id 



B„ 



eff 
9 



(7) 



with C = s%2f\VtbV t *\, and 



12tt 



d 
d 



(pB{X) 
1 — X 

<Pb{x) 
1 — X 



dx + 



dx 



>b(x 



X 



X 



dx 



dx I . 



After squaring the amplitude, and then performing the phase space integration over one 
of the two Dalitz variables, we get the differential decay width versus the photon energy E^ 1 



dT 2a 3 G| |T . ( . ;| . 2 



dE y 12 2 vr 4 



V tb V;f(C 2 1 +C 2 )(M Bq -2E J )E 1 



M Bq - 2E 

In numerical calculations, we use the following parameters [T5] : 



2M B d i 

+ °10 



(9) 



a = -±-, G F = 1.66 x 10- 5 GeV~ 2 ,^ = 0.4, f B = 0.19 GeV, 
"b. = 0.5, f Bs = 0.24 GeV, \V tb \ = 0.999, \V td \ = 0.007, \V ts \ = 0.041, 
M B o = 5.28 GeV, M Bs = 5.37 GeV, t b = 1.54 x 10" 12 s, t Bs = 1.46 x 10~ 12 s. (10) 



Table 1: Comparison of branching ratios with other model calculations 





Our Results 


Quark Model [7] 


light cone [H] 


BR(B S -> 7/^~) 


J--' +0.98 X iU 


4.6 x 10- 9 


1.9 x IO" 9 


BR(B S -> 7 e+e-) 


i.9;°f x io- 9 


6.2 x IO" 9 


2.35 x IO" 9 


BR(5° -> 7//+//-) 


0.65^;! x 10- 10 


6.2 x IO" 10 


1.2 x 10~ 10 


BR(£?° -> 7 e+e-) 


0.83;°;^ x io- 10 


8.2 x IO" 10 


1.5 x 10~ 10 



After integration of phase space, we get the decay branching ratios shown in Table [T] together 
with numbers from other models. The input parameters a, fs and CKM factors will only 
give an overall factor to the uncertainty of branching ratios, which can be obtained easily, 
thus we do not show them. The uncertainty shown in the Table [1] comes from the heavy 
meson wave function, by varying the parameter Ub = 0.4±0.1, and a>& a = 0.5±0.1. From this 
strong sensitivity, we know that the radiative decays can provide information or constraints 
on the B(B S ) meson wave functions. Surely more uncertainty can come from the next-to- 
leading order contribution which is difficult to estimate. From the table, we can see that 
our results are similar to the results from the light cone sum rule [8], but smaller than the 
constituent quark model [7J . If we interpret the second term of eq. ([HD as the inverse of the 
constituent quark mass, we will find that it corresponds to ma ~ 473MeV which is larger 
than that used in ref.p]. The simple picture of constituent quark model gives a larger result 
due to the smaller constituent quark mass. But the much smaller result for B° decays, which 
is even smaller than the light come sum rule one [8], is due to the smaller CKM factor \V t d\ 
used here. 

The differential decay width as a function of E 7 is displayed in FigJSJ Most of the 
contribution is with an energetic photon which is easier for the experiment. One may expect 
that through the mechanism of vector meson dominance [IE], long distance effects also 
contribute to the process B°(B S ) — > ■y£ + £~~ (£ = e, /x). Some detailed calculation jTTj indicates 
that: the long distance contributions are significant in the resonance region. It is probably 
not easy to derive short distance information from the dominant long distance contributions 
by total decay width. Therefore, the photon energy spectrum shown in Figj3]will be helpful 
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Figure 3: Differential decay rate for B S (B°) — > / ~f£ + £~ versus the photon energy E^. 
to distinguish those contributions. 

4 summary 

We calculate the rare decays B°(B S ) — > 7^ + £~ in SM. Utilizing the B (B s ) meson wave 
functions constrained by non-leptonic decays, the branching ratio is predicted at the order 
of 10 -9 for B s — > r yl + l~ and 10 -10 for B° — > / y£ + £~. There are possibilities to detect them in 
the LHC-b experiments, which could provide information on the B (B s ) meson wave function 
or new physics signal [T8] . 
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